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SUMMARY 


The problem of broadband noise generated by turbulence impinging on a downstream 
blade row is examined from a theoretical viewpoint. Equations are derived for sound power 
spectra in terms of 3 dimensional wavenumber spectra of the turbulence. Particular attention is 
given to issues of turbulence inhomogeneity associated with the near field of the rotor and 
variations through boundary layers. Both radial and blade-to-blade variations in turbulence 
statistics are treated rigorously without adding significant complexity to the final noise 
equations. Lean and sweep of the rotor or stator cascade are also handled rigorously with a full 
derivation of the relevant geometry and definitions of lean and sweep angles using 2 independent 
angle conventions. This report treats noise response of a cascade to turbulence with specified 
statistical properties; it is assumed that turbulence properties will be provided from CFD 
calculations or measurement. 

Derivation of the noise equations is made possible by the appearance of a flat plate, 
rectilinear cascade acoustic response theory by Professor Stewart Glegg of Florida Atlantic 
University. His theory includes a spanwise flow component (for treatment of the sweep effect) 
and a spanwise wavenumber component (in addition to the usual streamwise and gapwise 
wavenumber components in 2D theories). Glegg’s theory is for waves that are sinusoidal in 
space and time. This report provides the statistical analysis required to apply Glegg’s theory to 
the turbulence interaction problem. The final noise equations give the spectra of upstream and 
downstream sound power generated by turbulence interacting with a rotor or stator cascade. A 
specific recipe is given for dealing with inhomogeneity in terms of integrals of the turbulence 
covariance function (or of the turbulence spectrum function) at the cascade face. 

The general theory has been incorporated into the fan noise prediction system BFaNS. 
This requires, as a minimum, providing the circumferential and radial distribution of turbulence 
over the cascade inlet plane. However, use of the general theory is illustrated herein by 2 simple 
theoretical spectra for homogeneous turbulence. One is the Liepmann spectrum (for isotropic 
turbulence) and the other is an axi-symmetric turbulence spectrum from the literature. Limited 
comparisons are made with data from model fans designed by Pratt & Whitney, Allison, and 
Boeing. 

Parametric studies for stator noise are presented showing trends with Mach number, vane 
count, turbulence scale and intensity, lean, and sweep. Two conventions are presented to define 
lean and sweep. In the “cascade system” lean is a rotation out of its plane about the streamwise 
coordinate axis and sweep is a rotation of the airfoil in its plane. In the duct system lean is the 
leading edge angle viewing the fan from the front (along the fan axis) and sweep is the angle 
viewing the fan from the side (perpendicular to the axis). It is shown that the governing 
parameter is sweep in the plane of the airfoil (which reduces the chordwise component of Mach 
number). Lean (out of the plane of the airfoil) has little effect. Rotor noise predictions are 
compared with duct turbulence/rotor interaction noise data from Boeing and variations, including 
blade tip sweep and turbulence axial and transverse scales are explored. 

The report also serves as documentation of the equations in BBCascade and as a 
reference for the BFaNS code. BBCascade treats a single geometry and mean flow; BFaNS 
integrates the BBCascade equations over geometry and flow as they vary over radius of a fan. 
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SECTION 1 
INTRODUCTION 

Figure 1 represents the side view of a turbofan with an upstream rotor and downstream stator. 
Broadband noise is generated when turbulence interacts with leading and trailing edges of stator vanes 
and rotor blades. Three of the leading edge interactions can be analyzed with the theory ot this report 
and for each of these, the turbulence is intrinsically inhomogeneous or anisotropic, or both. For 
example, turbulence drawn into the inlet from the atmosphere or from flow over external objects 
interacts with the rotor. If the fan is stationary or moving slowly, the turbulent eddies are stretched in 
the flow direction so that the inflow turbulence tends to have an axi-symmetric form of anisotropy. In 
another mechanism, turbulence in the boundary layer of the duct interacts with the rotor. Since this is 



Figure 1 Noise generation in a turbofan by Fiqure 2 Inhomogeneity of rotor wake at entrance 

turbulence in rotor wakes impinging on stator. to s tator.' Test data from reference 1 . 

confined to a thin layer, it is highly inhomogeneous in the radial direction. Finally, wakes from the 
rotor blades impinge on the stator. The turbulence field behind a rotor is shown in Figure 2. It also is 
highly inhomogeneous with turbulence statistics varying periodically in the gapwise direction. A 
fundamental prediction theory is needed for these complex mechanisms that is sensitive to turbulence 
properties (scales and intensities) and that rigorously addresses the issues of inhomogeneity and 
anisotropy. Furthermore, with the current interest in lean and sweep of stator vanes and rotor blades, 
the theory should treat those geometries rigorously, also. 

In this report, we develop theory for the spectra of upstream and downstream sound power from a 
rotor or stator cascade as an integral over the wavenumbers of the incoming turbulence field. The only 
simplification regarding representation of the turbulence is use of Taylor’s hypothesis (i.e. the frozen 
gust assumption). This reduces the turbulence representation from a 4 dimensional wavenumber- 
frequency spectrum to a 3D wavenumber spectrum and results in manageable noise equations. The 
principal simplification regarding geometry is use of a rectilinear, flat plate cascade model. At low 
frequencies, this may be problematical; however, at mid and high frequencies where noise generation 
is dominated by small scale eddies, the rectilinear cascade approximation should be acceptable. 
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Some other aspects of this work have been published previously. Reference 2 examined the 
problem of measuring turbulent inflow to a stator for the purpose of making noise predictions. Two 
probe measurements are required and guidance for the required data processing and probe separation 
were given. Reference 3 showed in detail how the turbulence spectrum and the acoustic response 
spectrum work together in a stator noise calculation. Parametric studies showing trends with Mach 
number, vane count, turbulence scale, etc. were also given. Reference 4 is a greatly reduced version of 
this report that adds effects of lean and sweep to the previously published results. 

Tins work has „ cn enabled by the appearance of a harmonic cascade theory by Glegg (ref. 5) that 
this report adapts to the broadband problem. Glegg, too, has published broadband applications using 
his cascade theory. Reference 6 is a recent example. The main additions of the present report arc the 
focus on lean and sweep plus treatment of turbulence inhomogeneity and anisotropy. It also serves as 
documentation for the code BBCascade, which treats rotors and stators in a unified analysis. 

Analysis herein is confined to noise from a cascade with constant geometry and flow properties in 
the spanwise (or radial) direction. It is a useful analysis tool as is; however, for fan noise diagnosis 
and design it has been incorporated in Morin’s BFaNS code (ref. 7). In BFaNS, geometry and flow for 
the leading edge sources on the rotor and stator are integrated over the fan radius using the theoretical 
equations of this report. 

In the present report; Section 2 is a review of Glegg’s harmonic cascade theory. Section 3 
outlines the geometry of lean and sweep and gives the equations for transforming coordinate systems 
and wavenumbers from duct coordinates to Glegg’s cascade coordinates. Section 4 is the derivation of 
the broadband noise theory. Section 5 presents the 2 turbulence spectra currently in use: one for 
isotropic turbulence and the other for axi-symmetric turbulence. Section 6 non-dimensionalizes the 
working equations and provides documentation for code BBCascade. Section 7 shows how the noise 
prediction equations work with illustrations of turbulence spectra and cascade response spectra. 
Section 8 shows the capability of the broadband theory via comparisons with test data from Pratt & 
Whitney, Allison, and Boeing test fans and by parametric variations to explore effects of turbulence 
scale, flow Mach number, stagger, solidity, vane lean & sweep, vane and blade count, and rotor sweep. 
Appendices A and B give details of coordinate and wavenumber transformations required for the 
analysis. Appendix C provides a list of notation. 
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SECTION 2 

BACKGROUND ON GLEGG’S HARMONIC CASCADE THEORY 


Previous noise models for turbulence/stator interaction have been based on simplifying 
assumptions some of which are no longer necessary. Some examples: treatment of the blades or vanes 
as isolated airfoils (as opposed to cascades), 2 dimensional theory, incompressible flow for airfoil 
loading response, and compact (as opposed to distributed) radiation theory for the acoustic sources. A 
recent analysis by Cilcgg' re ' overcomes these limitations but still is tractable enough for calculations 
on ordinary computers. Furthermore, it can be adapted to leaned and/or swept cascades by virtue of a 
spanwise velocity component included with the mean flow. 

As shown in Figure 3, geometry is constant in the z direction and the background flow is 
uniform: U =(U,0,W). Airfoils are unloaded flat plates. Cascade gap, chord, and stagger angle are s, 

c, and x • 



The unsteady flow is harmonic in space and time with upwash given by Equation 1 

/ i(y n x+ay+vz-a>t) 

w(\,t) = w 0 e y '° ' 


0 ) 


This represents a plane wave that is harmonic in time with frequency a> / In and upwash complex 
amplitude w Q . It is also harmonic in space with x, y, and z wavenumbers equal to y 0 , ct, and v. By 
use of Wiener-Hopf analysis, Glegg derived an equation equivalent to the following for the velocity 
potential of acoustic waves scattered by the cascade in response to the input wave of Equation 1 

2 y 1 c'l-A (x-yd/h)+((T-2nk)y/h+vz] p -icbt (2) 


+ 7tW n C 

<fr(x,t) = ± 0 


P<e 


where the +/- signs refer to upstream/downstream going waves and 
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N 

11 

£k = -inlr 


*Ik =~fk sin Ze ±c 0S Xei 

/*•«-/< 


= 03 ~ IV V 

d / fih = kM + 77 ^ 

(W/?) 2 k = (o g / (aj3 2 ) 

f\ -{<7- Ink + tcMd) / s e 


( 4 ) 

( 5 ) 

( 6 ) 


and a = 7 0 d + ah is the interblade phase angle. D is the Fourier transform of the discontinuity in 
potential across the blade and wakes (in the form of an infinite product) and is the major result of 
Glegg’s derivation. The velocity, pressure, and density perturbations associated with the acoustic wave 
can be obtained from Equation 2 via 

u = V 0 p = - p 0 D(f> / Dl p’ = p / a 2 (7) 

where p 0 and a are the ambient density and speed of sound. We adhere to Glegg’s notation closely. 

However, he used a>' for the radian frequency (our oj ), a) for the shifted frequency (our ojg), and 
m for the scattering index (our k). Also, the D function of this report (and of Glegg’s code) is non- 
dimensional; to obtain the D function of Glegg’s report, multiply the non-dimensional version by 

w„C. 


The formulation above gives the acoustic waves scattered by a cascade for a single planar wave 
input per Equation 1 . Scattering index k runs over an infinite range but, as usual in this kind of 
formulation, only a finite number of waves are cut on (propagate undiminished); the remaining waves 
decay exponentially and, thus, carry no acoustic energy. Cuton is governed by the argument of the 

I 2 2 

square root -\jK e - f k ; when the frequency is high enough, the argument is positive and the waves are 

cut on. The exponential dependence of the space and time variables in Equation 1 permits treatment of 
any inflow field via standard Fourier transform methods, as shown in Section 4. 
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SECTION 3 


COORDINATE TRANSFORMATIONS AND 
GEOMETRY FOR LEAN AND SWEEP 

This report requires several coordinate systems for the different parts of the analysis. One is 
aligned with the fan duct for the lean/sweep definition, one is aligned with the mean flow and the fan 
radius to define a turbulence symmetry axis, and one is aligned with the cascade geometry to apply 
G1 egg’s theory. Furthermore, rotating coordinates are needed for the rotor and stationary coordinates 
for the stator. The required transformations are treated in detail in Appendix A. In this section, we 
define the transformations and lean/sweep conventions and we summarize the most important 
geometrical relations between chord, gap, and stagger in the duct system and in the cascade system. 
We also define lean and sweep in 2 different conventions. 

The coodinate systems connecting duct coordinates to Glegg’s coordinates are shown in Figure 
4. Stagger, lean, and sweep angles are defined via an ordered sequence ot rotations about coordinate 
axes as follows. The sketch at the top represents the stator or rotor cut at constant radius and 
unwrapped onto the plane of the paper. The z d and z 0 axes are coincident with each other and with a 
fan radius. Mean flow U Q is aligned with the direction and has no radial component. The 
following explanation of Figure 4 is in the context oi vanes in a stator but it applies to rotors as well. 
In the transformations, the reference vane (vane #0, which passes through the origin) and vane #1 are 
handled differently. To visualize the rotations, think of rotating the reference vane and then rotating 
the coordinate system to re-align the z axis with the leading edge of the reference vane. For the other 
vanes, think of them as having hooks on their leading edges attaching them to the yd axis with a 
separation g that does not change. Thus, throughout the lean and sweep rotations, the vanes always 
penetrate the constant radius plane (z 0 = 0) at the location and orientation given by the gap g and the 
stagger 0 at the top right of Figure 4. 

Lean is an out of plane of the vane rotation about the x 0 axis through angle y/i . Note that lean 
does not change the component of inflow mean velocity normal to the leading edge. Sweep is then an 
in plane rotation about the y ' axis through angle y / s . This does reduce the flow component normal to 
the leading edge by the factor cos y/ s . Thus, transformation from duct coordinates {x d , yd , z<d t0 
Glegg’s cascade coordinates (x, y, z) requires an ordered sequence of rotations (stagger, lean, sweep) 
about the z d , and y ' axes. The net transformation, derived in Appendix A, is represented by 


X 


x d 

y 

= [2]x 

yd 

V 


fd. 


where 



cos 6 cos y/ s + smOsmy// sin if/ s 
-sin 6 cosy / 1 

cost?sin y/ s - sin^sin y/ / cosy/ s 


sin Ocos y/ s - cosflsin y / / sin y/ s 
cos 0 cos y / 1 

sin 0 sin y/ s + cos 9 sin y/ [ cos y/ s 


-cosy / 1 sin y/ s 
-sin y / 1 
cos y/j cosy/ s 


(9) 
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STAGGER ANCLE 0 
Vane sections at constant radius 
unwrapped onto jq/, ya plane 
(duct coordinates). 
x 0 , Vo coordinates aligned with 
vanes at stagger angle 0. zj and 
z 0 axes aligned with leading 
edge of referenee vane. 


Other views in the x 0 , y 0 , z 0 
system. For the lean and sweep 
operations to follow, consider 
the vanes attached to the ya axis 
with “hooks'”. Coordinate 
systems follow the reference 
vane; vane #1 stays parallel to 
ref. vane (hooked at y ( i = g). 


V 0 


y 0 


__X() 


LEAN ANGLE y/, 

Lean the reference vane about 
the x a axis through angle y/\. z 
axis follows leading edge, x 
axis coincides with x a axis. 
Because of the “hook”, the gap 
in the xj, yj system is not 
changed. 


y 0 


SWEEP ANGLE y/ s 

Sweep the vane in its own plane 
through angle y/ x . Velocity 
vector does not change, c,/ is 
the chord measured at constant 
radius in the fan. 
c - Cd cos y/ s 


V 0 


The x, y, z coordinates are now 
Glegg’s cascade coordinate 
system. 







Figure 4. Sequence of rotations to relate cascade coordiates (x,y,z) to duct coordinates ( Xd,y d ,Zd ) and 
intermediate systems. 
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Relationships between the velocity components in cascade-fixed coordinates are easy to deduce 
from the figure 


U = U 0 cosy/ s 
W = U 0 sin y/ s 


where U 0 is the total mean inflow velocity (in cascade-fixed coordinates). 

If we say that the geometry descriptors in the duct system are chord cj, stagger 0 , gap g , lean 
angle (/// , and sweep angle y/ s and in the cascade system are chord c , gap s , stagger angle % > ancj 
sweep angle y/ s = axc\m(W/U), then the connection between the 2 systems is given by the following 
equations. The relation between the chords c and cy is simply 

c = Crf cos y/ s ( * ’ ) 

where cj is chord measured at constant radius (along the flow direction). Appendix A shows that 
Glegg’s stagger dimensions in Figure 3 are given by 

d = 0\2 8 = ( sin $cos - cos 0 sin y/j sin y/ s )g ^ 2 ) 

h = Q 2 2 g = ( cos ^ cos ¥1 )£ 

From these, his gap and stagger angle are 

5 = 4cfi+h 2 = g^Jo^+Qh 


tan* = 




(14) 


We can also define lean and sweep angles via a convention more suitable to viewing from the 
point of view of mechanical design. Consider the sketch below where the cascade is viewed along the 
fan axis for lean (on the left) and from the side for sweep (on the right). 



Figure 5. Lean and sweep angles defined via front and side views of stator 
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Here the lean and sweep angles have hats to distinguish them from the angles used in Figure 4. The 
relationship between the 2 systems is derived in Appendix A and is repeated below. 


tan i // 1 = cos 6 tan y j - sin 0 tan y s 

tan y/ s ~ cos y / (sin 0 tan y t + cos 0 tan y s ) 

Equations 1 3 must be solved sequentially. Thus, given lean and sweep angles in the “duct system” 
above, they can be converted to the “no hat” system for use with Equations 9-12. The inverse of 
Equations 15 giving y j and y s as functions of yi and y s can also be found in Appendix A. 
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SECTION 4 


DERIVATION OF EQUATIONS FOR SOUND POWER SPECTRUM 

In developing the theory of this report, the derivation was made as general as possible. The result is a 
single analysis that applies to noise from rotors and stators with lean and sweep. Turbulence is 
represented by a 3D wavenumber spectrum defined in a manner to include inhomogeneity and 
anisotropy. It will be seen that the general expression for acoustic intensity is too cumbersome for 
numerical work. Intensity varies with radial and gapwise position of the field point in the duct. 
However, when intensity is integrated (or averaged) over the duct area to compute sound power , the 
formulas simplify dramatically and tractable equations are found that treat inhomogeneity rigorously. 
It turns out that averaging in the field coordinates permits an averaging over the turbulence properties 
in the gapwise direction and an integration of the turbulence in the radial direction. Thus, the 
turbulence spectrum needed for input to the theory must represent the gapwise average of the field; for 
stator inflow, this is what would be measured by probes at the stator inlet with time averages over 
many rotor revolutions automatically giving the correct gapwise averaging of the turbulence field. 

The sketches below show the general case analyzed and how it can be applied to rotor and stator 
problems. Rotors can be analyzed with axial inflow or swirling outflow, but not both. 


M To f St at M r 



Figure 6. Mean flow viewed in constant radius “plane". General case with swirl and rotor rotation plus 2 
specialized applications: stator with swirling inflow and rotor with axial inflow. 


4.1 Generalization for Random Inflow 

Glegg’s formulation described in Section 2 was written for 3D planar, harmonic waves. His theory can 
be extended to an input upwash w(x,t ) with any waveform and spatial distribution via the Fourier 

transform 

w(x,0 = jj ^(K,5)e* K '*~ a) dK«/S (16) 

where here, and throughout this paper, integration limits are from -oo to oo unless otherwise specified. 
We placed a tilde (~) on the co to denote frequency in the cascade frame. Most of the derivation 
below uses this source frequency; at the end this is shifted to observer frequency in the stator frame to 
obtain the desired acoustic spectrum. Vector wavenumber K is shorthand for (y 0 ,a,v) and x 
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denotes position in a coordinate system fixed to the cascade being analyzed (rotor or stator) and 
replaces Glegg’s coordinates. For later use, the inverse of Equation 16 is 


JP(K,<y) = w(x, t)e ,(K x M) dxdt (17) 

With the application of Equation 17 to Equation 2, we can generalize Glegg’s potential immediately to 


?*(!,/) = ±-^- 


JjM'fK.S) £ 


V Ah -Ak)_ ‘f l + -?~ L )y+v-< rf] 

! 2 " -2 

A-.-OC ^fC e - j k 


dKdco (18) 


where the tilde on means that the expression for potential applies in the cascade-fixed frame. 


Tra nsformation to Stationary Coordinate System 

We are interested in computing the spectrum of sound power scattered by the cascade. The starting 
point is the expression for the acoustic energy flux vector applicable to waves in a uniformly moving 

medium with mean properties given by the density p 0 , speed of sound a, and velocity U (ref. 8) 

/ \ 

1= — +U u {p a u + \}p’) (19) 

\Po J 

This is the time dependent power per unit area in terms of the acoustic pressure, density, and velocity, 
p, fd , and u. Since it depends on the mean velocity U , it must be computed in the stationary 
reference frame (fixed with respect to the duct and observer) whereas Glegg’s expressions for the 
scattered acoustic field apply in a frame fixed to the cascade. The coordinate shift, expressed in duct 
coordinates (with subscript d in Figure 4), is y d = y d -Q Rt . Parameter R is the effective radius of 
the source and is also used for scaling later in the derivation. The transformation matrix Q defined in 
Section 3 rotates this to Glegg’s system so that the coordinate increments for the transformation, 
5 = x + Ax , to the stator frame are 


1 [ 0 1 \Qn 

Ay = Q x -Q R( = - Q 21 QR t ( 20) 

A*. - 0 .032. 

Adding these deltas in the exponent of Equation 1 8 results in the frequency shift 

03 k = ~[(~0I2 + ^Q22) Xr k + 022(^~) + 032fP^ (21) 

However, by use of Equation 12, the first 2 terms cancel, leaving 

<°k = ~[(^p) + v]Q* ( 22) 

An alternative expression for ct)k can be obtained by inserting (r=y 0 d+ ah for interblade phase angle 

and by noting that 2nR/g = B is the blade count in the rotor (or vane count in the stator) and by again 
using Equation 12. 

<°k = -[01 27 o + 022« + 032 + kBQ ( 23) 
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In the stationary frame, the acoustic potential is now 


fi ± =± 


TCC 




fM K - 


CO 


>1 

A~-qo 


ciix 4 ) 

*■« - /* 


i|-4*+Ul )'] (/K ( 24) 


expressed as an integral over the turbulence wavenumber vector K. and the turbulence frequency (o . 
Recall that a> is turbulence frequency in the rotor frame-, this form is convenient for now but will 
eventually be replaced by one with integration over observer frequency to obtain the desired spectrum. 
Coordinates x,y, and z are the same as Glegg’s for stator problems and are parallel to Glegg’s (but 
moving) for rotor problems. 


4.3 Acoustic Intensity 

The expression for power flux in Equation 19 can be written in terms of potential using Equations 7 as 


follows 


I = -Po 


dfi 

dt 


V fi- 


ll Dfi\ 
a 2 Dt) 


(25) 


For convenience in the ensuing manipulations, we can rewrite this as 


-Po 


dfi 


dt 


Vfi- 


U Dfi 
a 2 Dt 


(26) 


where * denotes complex conjugate. This is permitted because the potential expressed in Equation 24 
is pure real as a result of using a double sided summation and double sided integrals. (I he imaginary 
parts in the upper and lower halves of the summation and integration ranges cancel.) There is no need 
to imply “real part of’ with the notation. Thus, since fi is pure real, it follows that fi= fi . 


To find the acoustic power leaving the cascade, we will compute the component of intensity normal to 
the cascade 

/* = n* I (27) 

and then later integrate it over the frontal area of the blade row. Here, n* are upstream/downstream 
unit vectors normal to planes containing the cascade leading edges or trailing edges 


+ 

n” 



(28) 


Formally, Equations 26 and 27 must be evaluated in the coordinate system fixed relative to the 
observer (or stator) and it is helpful to establish the role of blade row motion in the required dot 
products and derivatives. For example, dfi/dt yields — which explicitly contains the 

frequency shift due to the relative motion of the rotor. n - V is the same in the rotor and stator 

frames, n 1 • U can be evaluated in the rotor frame since the relative motion is perpendicular to n . 
Finally, since Dfi/Dt is the convective derivative, it must give the same result regardless of coordinate 

system used for its evaluation. Hence, n 1 ■ ( V0 — tTg) as we ^ evaluated in the rotor 

a~ Dt 

frame since the process is easier. The result is 
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T ± ± Po * C 


J j | W *(K< a>W(K', S')£ X (5 - ) 


(4)^(4) 


Ke-fH*'e-fk 


x {'- A? + f (A? j + **=%—) -7 ^ (S' + UA? - V'W)} ( 29 ) 

xe , '{-(4 -■*£ M+[(4 i i + ~h~ >b'+( Hw’-o4. )]/} 

xdK dK'dcodco' 

where now the sums and integrals are over the limited ranges of k and k! and wavenumbers 
corresponding to cut on (energy carrying) modes. Glegg’s found that the term in curly brackets { } is 

± Bs , I 2 ~ 

equal to - f k > , which simplifies Equation 29 considerably. Since we are dealing with 

.v 

turbulence statistics, we take the expected value of both sides, bringing Equation 29 to 


2 4 

j±\ = Po£ g 


— jjjj( »’’(K,5)»(K',S') }££ 

AA<J A A' 


cirtD\4)w $) 

-./ft" 


-i{ -(a* ->1*. )*+[U* i+—fr )-i A k' <l+~-j,- )}y+(v-v')~-\{a)-6> k )-{w'-(o' k . )j/j 


xdK dK'doidd)' 


Equation 30 represents power per unit area of cascade face. It is still a function of the space variables 
x v, z and time. For further simplification, we will analyze the upwash spectrum 

(w (K,ry)JF(K',<y')^) in the next section where we show that the expected value operation produces 

delta functions that enable some of the integrals in Equation 30. 

4.4 Reduce Upwash Spectrum to Standard Form for Turbulence Spectrum 

In this section we examine the flows typical of fans and develop an expression for 

(w (K.,d))W(K' ,a>')) in terms of a 3D turbulence spectrum while retaining as much generality as 

possible with respect to inhomogeneity and anisotropy. We will address the most difficult type of 
inhomogeneity, namely the rotor wake entering the stator, and the other applications (duct boundary 
layer/rotor interaction, etc.) will follow. Figure 7 (at the left) is a view looking radially in the fan 
indicating the discrete nature of the wakes and the stator vanes. At the right is an axial view indicating 
a band of turbulence, so that we must deal w ith inhomogeneity in the radial direction. 



Figure 7. Sketches showing types of inhomogeneity in wakes of rotor entering the stator. Left: view along 
radius. Right: view along fan axis. 
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The procedure for evaluating < IV* W) is to substitute the inverse transforms from liquation 17 and 
perform some of the integrals using generalized function theory. Thus, 




At this point we adopt Taylor’s hypothesis, treating the turbulence as a frozen gust pattern convected 
past the cascade. We use the notation w for the frozen turbulence field, which is time-independent in 
the fluid-fixed frame. Its convection in the x 0 direction (see Figure 4) would normally be expressed as 
w(x,t)= w(x 0 -U 0 t,y 0 ,z 0 ). However, at this point we employ a device for the treatment of 
inhomogeneity. We insert a random time offset T representing the fact that we do not know what 
point in the wake passing cycle corresponds to t — 0. Thus, we write (noting that z 0 zj) 

w(xJ) = w[x 0 -U 0 (t-T),y 0 ,z d ] (32) 


Random variable T eventually disappears from the equations. Now we display the coordinates y 0 , 
Zd explicitly in the integrals 


W*W 


JJJiJJJJ ^ vT ^ X ° “ U °^ ~ r) ’ y ° ,Zd ^ X '° ~ U ° (r “ r) ’ y °' z ' d ^ 


(2/r) ? 


(33) 


To proceed, we first shift back to fluid-fixed coordinates via x Q = x + U 0 t and x' 0 = x' + U 0 t ' . 
Then we change from using 2 independent points to represent the function (e.g. y n and y’ 0 , etc.) to an 
“anchor” point and a spacing (e.g. y 0 and s y ) via the following coordinate shifts. 

X ' = X + s x 

y'o = y 0 + *y ( 34 ) 

z' d =z d + s z 

t' - t + T 


In Equation 33, k x ,k y , k z are associated with the x 0 ,y 0 , z d axes. With these substitutions, 

w(x +U 0 T,y 0 ,z d )w(x +s x +U 0 T,y 0 +s y ,z d +s z )} 
xe -i(k’ x s x +kySy+k’ z s z ) e i(k z -k' z )z d e i{k x -k' x )x e ‘(k y -k’ y )y 0 

xe ‘[(k x U o -ci))-(k x U 0 -a> )]t e ~i(k x U 0 -a) )t ^ dz d ds x ds y ds z dt dt 

There follows now a critical step regarding inhomogeneity of the turbulence. The expected value 
of the velocity product on the right side of Equation 35 is expressed as a function of the fluid-fixed 
position variables v and y 0 , as if we knew the / = 0 point for the velocity pattern. However, recall 
that we inserted the random variable T to represent the fact that we do not know this point. Hence, if 
the expected value symbol ( ) denotes an ensemble average, then the ensemble is all values of T in a 
blade passing period. Since these are equally probable, the average is independent of x . Once we 
have averaged out this jc dependence, it is easy to see that the expected value (tw) is also 
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V ' 

TURBULENT 

Vi / ■ A -' 

«** ■ V l -' '1/ v 

TURBULENT 

V'-- --\-i 


Figure 8 Sketch showing degenerate case not 
covered by analysis 


(wiv) as the velocity covariance function 


independent of y a (but still dependent on the radial 
coordinate z d ). Note that, if the 2 blade rows were 
not in relative motion, then the left side of Figure 7 
would look like the sketch in Figure 8. In this case, 
the ensemble average would not be independent of 
the transverse coordinate and the ensuing analysis 
would have to be modified. 

In light of the above discussion, we recognize 


-H(s, Z d ) = (m'(x)w'(x + s)) 


(36) 


where the overbar on 9? denotes the averaging just discussed. It is explicitly shown to be dependent 
on Zd and, by implication, independent of x and y n . With this notation 


W 




-i(k' x s x +k's -t k' 2 s. ) 


ds x ds y ds z 


e U z -K)Zd dzd 


(2 n)- 


-jeM'-kxV&ji 




dy„ j e ' 1 


i[(k x U 0 -cb)-(k' x U o -a)')]t j ( -i(k' x U 0 -(b')T 


37) 


dr 


The expression in square brackets | ] is an extension for radial-dependence of the standard 3D Fourier 
transform (Ref. 9) defining the turbulence spectrum function 

<j> 22 (K',z d ) = JJj9i ( s,z^)e- ,K ' V S ( 38) 

The integrals on the bottom line of Equation 37 are all delta functions so that we arrive at 

(w w ) = Y n 8(k x - k’ x )6{k y - ky)S(a> - d)')S(k x U 0 - d)’)j ^ 22 (K',z^)e' (A ‘ )Z(t dz d (39) 

Because the x , y 0 ,z d and x d ,y d ,z d systems share an axis, S(k x - k' x )S(k v - k' y ) specifies a point in 

the k x , k y plane that can just as well be denoted by 5{k xd - k' xd )5(k yd -k' yd ). Hence, Equation 39 
can be written alternatively 

{ W ' W ) = &( k xd - Kd)S{k yd -k' yd )S{a-a')5(k' x U 0 -a')]^ 22 {K\z d )e i ^- k ^dz d ( 40) 

To simplify notation, we have left the argument of the last delta function unchanged but note that 

k' x = k xd cos 6 + k' yd sin 6 (41) 
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4.5 Integral of Intensity for Sound Power 

We can now insert the expression, just derived, for the generalized upwash spectrum {W W) into 
Equation 30 for intensity and then integrate over planes parallel to the face of the cascade according to 
n ± =|(/ i ) dA to get sound power upstream and downstream. 


Combining Equations 30 and 40 and performing the co' integral give 


2 A 


j i «**/ - *;</ wv - k 'ydWAo - ®)J te(K’. zjV ik -- k ' ! * 


X ZE ( G>-G>k ) 

A A' 


(5^d‘(4)D(4) 




(42) 


xc 


-d-(4 - a? )*+l(4 f . „ 


±' | a' -lnk* 


')b- +( v-v')z +(®A-«E)O i/Kt/KW ~ 


Integration over the face of the cascade is best done in the “hat” 
coordinates of Figure 9 in which the j>,£ plane contains the leading 
edges of the cascade and y is the tangential coordinate (at constant 
radius) coincident with the w axis. Appendix B shows that rotation 
of the x, y, z system into the “hat” system changes the long 
exponential in Equation 42 to the following (in which we have used 

the fact that k y j = k vc j ) 



-i{(* i -*5;. )- 

e 8 ” 


(43) 


It will be seen below that the form of the axial wavenumber of the response waves k x j ( need not be 
known for the derivation; we only need know that it depends on k and the excitation wavenumber. 
Recognizing that the primed wavenumber vector differential can be expanded as 
dK' - dk' xc jdk' y( jdk' zc j , we can replace the long exponential in Equation 42 with Equation 43 and 

perform the k’ yci integration via the delta function 8(k y d - k ' y J) 



.2 A 


1 1 S(k >d - VdWWo - S >1 


k k' 




■vd-fk 


2 K 


A{{k? k -k a xk .)i-— (k-nM(kzd~hd)--^(k-kA]zHk-k')BQt} 


xe 


dK dk' x( jdk' zc jd(o 


( 44 ) 
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Now, we can do the first step in the area integration, namely the integration over the j> coordinate. 
To obtaine the desired sound power, we integrate in this direction over B blades (or vanes), or for a 
distance of Bg where g is the gap between blades. This integral involves 

{k-k')y 

dy=Bg5 kk . (45) 

0 

Kroncker delta 8 k k' enables the k' summation, leaving 


Bg -i 2 * 


!('■>*' ■ J S£r mi 8(k xd - k' xd )5{k' x U 0 - S)J M K\z d )e i(k :- k * )Zd ciz d 

± ; ± * ± (± a a (46) 

x]T (S - CO k ) g ~ ,T( *** -Kk )k+(k 2 d -k' :d )i\ dK dk'^dK.idd) 


*xd aK zd c 


< ~ fk 


Note that the time dependence has disappeared. We want to do the k' xd integration next via the delta 
function S(k xd - k' xd ). However, the wavenumbers in the exponential are in a different coordinate 
system. T he rotation for the wavenumbers is the same as for the coordinates shown in Figure 9 from 
which we determine 


k-zd k^d (k x d k xd ) sin y/ s + ( k zd k~ d )cosi// s (47) 

since k yd = k vd is common to both systems. Now we can do the k' xd integration with the result 



n 2 4 

Po ^gn c 
2 nfiss e 


111 


-S)|fe(K ’,z d y {k ‘- k 'o--«dz d 


X 


I 


k 


(do-co k dk'jdTo 

i K ~e ~ fk 


(48) 


And we complete the area integration of intensity to get sound power via 


n ± 


= J[J (!*)$]& 


(40) 


We interpret the outer integral as running from -qo to +oo so as to capture all of the radiated energy. 
This integral involves only an exponential and thus produces another delta function 



cos yt s 


(50) 


The term cos (j/ s in the denominator is in the direction of sweep increasing noise. It appears here 
because sweep increases the length of leading edge exposed to a band of turbulent flow. Later, it will 
be seen that this undesired tendency will be overcome by the reduction in effective Mach number due 
to sweep. Equation 50 enables the k \ d integration leading to 
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(51) 


J-J± _ Po Bg^ J- 


\ _ \Cl fi('k'k) 

-f f S(k x U 0 - <y)j <t>22^- z d)^ z dZ_j (<»-<«/() p^-=^ 

.v A- 4 K e~Jk 


f:!ss e cos y/ 


dKdd) 


Note that the exponential has disappeared entirely, including the term (k xk ~ ^xk' )* • 1 * 1 ' s occurs 

because we have arrived at the point where K=K' and k—k' so that k xk = k xk ' • 

Equation 51 has gap g in the duct system and gap .v in the cascade system. This can be 
eliminated by using the result from Appendix A which shows that 

8 I (52) 

.s’cos y/ s cos y/ / cos y/ s 


so that, with a change in order of integration, 


U ± = gp ^ - C - f f \MK,z d )ik d J^S(k x V 0 -at){a-a> k y-^ 

P$e cos yf i cos y/ s k ^ 


I 


.0 


U do) dK 


(53) 


■/* 


The integration order can be shifted without concern for integration limits if we take all the ranges to 
be -oo to +oo and consider the integrand to be zero whenever the square root in the denominator 
becomes imaginary (i.e. for cut off waves). Equation 53 represents the upstream and downstream 
sound power integrated over all frequencies rather than the sound power spectrum. Recall that the 
frequency integration variable 5 is source frequency rather than observer frequency. This can be 
shifted to an integral over observer frequency co via 

d)-Q) + co k (^) 

which is now permitted with the a> integral inside the K integral. But, then, with the infinite limits, 
the integration order can be changed back again so that the power is in the form 


IT 


-J. 


+00 > 
n- a dco 
—00 


where 


n£ = 


ki d ( 4) 

„ j [J ^22( K ' z d frfcrf ]Z 5 (k x U 0 1 ~ ~ 2 

/3s e cos y/ j cos y/ s ^ LJ } k ^ e ~Jk 


Bp 0 K 2 c^a> 


dK 


(55) 


(56) 


Recall that the Doppler frequency shift from Equation 22 is 

®*=-i( 2: f i )+e32''P* < 57) 

Equation 56 is the spectrum of sound power distributed over observer frequency, as desired. It is 
expressed as an integral over the wavenumber vector K of the turbulence. The final delta function 
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still needs to be dealt with to reduce the wavenumber integration from 3 wavenumber components to 2, 
but that step is postponed. 

Our current form for the spectrum is not efficient for computation because Glegg’s D function is 
inside the integration and summation. Another shill in integration variable can be used to move D 
out from under the sum at the expense of moving the turbulence spectrum under the sum. This is 
advantageous because <f > 2 2 w '^ normally be computed from a simple algebraic expression whereas D 
requires considerable computational effort. Note that dK can be expressed in various wavenumber 
coordinate systems but, for the present purposes, it is considered to be dK^d^ da d v. Also, note in all 
ol the expressions for the acoustic response functions and in the frequency shift (o^ that a always 
appears with the scattering index k in the combination ah - Ink . Thus, if we shift the a 
wavenumber variable to a' via 


ah — 2nk = a'h 


(58) 


we achieve the desired result 


ps e cos if / 1 cos y/ s : J J 




Pe~.fo A 


X [Jfc 2 (K ,zj)dz ct ^ dy 0 da’dv (59) 


where the scattering index k has disappeared from the freqency shift: 

°>o - ~(Q\ 2 V o + Qn a ' + 032 v)oy? 


(60) 


The delta function cannot be used for integration until the wavenumber k x is expressed in the /„ a' v 
system. Appendix B shows that 


k x ~ Y o cos^ v + vsin y/ s (61) 

Substitution of this and Equation 60 into the argument of the delta function gives, with some re- 
arrangement 


5 -> 5[(U + Q ]2 QR)y 0 ~(co-Wv) + (Q 22 a' + Q 32 v) QR] 
This is functionally equivalent to 


(62) 


1 


u + q ]2 qr 


s 


y 0 


( co - Wv)-(Q22a' + 032 
U + 0j2 kiR 


(63) 


This enables the y 0 integration, so that finally 


n* = 


Bp 0 n 2 c A co 


II 


k loud 


fis e (U + Q ]2 QR)cosy// cosi// s ^2 _ J 2 ^ 


X [jfe (K,z d )dz d da' dv (64) 
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and, in evaluation of Glegg’s acoustic response function, the values used for y„ and (b are 


( co — Wv) — (Q 22 a> + G 32 

r °~ Z7+ QuQR 

and 

(b = CD-(Q\2Y o + Qll a ' + V )^R 


(65) 

( 66 ) 


Finally, recall that in evaluation of the turbulence spectrum that the a wavenumber must be shitted 
according to 


a~a' + 2nk/h (67) 

The equations just derived provide the transfer function from the turbulence spectrum as input to the 
output spectrum of sound power. In the following sections, we establish forms for the turbulence 
spectrum and non-dimensionalize the equations for coding. 


NASA/CR— 200 1 -210762 


19 



. 






SECTION 5 

TURBULENCE SPECTRA FOR BBCASCADE 

The theoretical acoustic spectra derived above, when applied to inhomogeneous, anisotropic 
turbulence, require that the turbulence spectrum be averaged properly in the gapwise direction. 
If this kind of test data is available, then it should be used. However, for predictions based on 
theory, turbulence spectra are not yet available for this general kind of flow. Out of necessity 
BBCascade uses spectra for homogeneous turbulence available from the literature. Thus, the 
^ 27 of Section 4 will be replaced simply by <h.2 , as given in this section. Two options are 
provided: one for isotropic turbulence and one for axi-symmetric turbulence. For the isotropic 
case, we use Leipmann’s spectrum (ref. 10). For axi-symmetric turbulence, we use theory given 
in a paper by Kerschen and Gliebe (ref. 1 1). This section presents the algebraic forms ot the 
spectra, shows how to rotate wavenumber coordinates to make the functional dependence of the 
spectra compatible with the noise equations, and normalizes the spectra for use in the computer 
code. Both spectra are coded in normalized form <&22 based on the local radius R and the 
mean flow velocity in the stationary frame U sta , . The normalization is defined by 

022 = ^^22 < 68 > 


Wavenumbers will also be normalized based on R so that K/—k/R , etc. 


5.1 Liepmann Spectrum 

This is one of the standard 3D spectra often used for acoustic and unsteady 

. .. 2m^A 5 *i +*3 

22 n 2 [l + A 2 (fc, 2 +*2 +fc 3 )I 3 


air loading analysis 

(69) 


2 

It applies the velocity component in the “2” direction whose mean square value is w . A is 
the integral scale of the turbulence. Because this is a spectrum for isotropic turbulence, it does 
not matter how the coordinates are oriented in space except that the “2” direction must be 
normal to the airfoil surfaces (Glegg’s y direction). Since we are free to select the orientation of 
the other axes, we make the most convenient choice and align the turbulence coordinates with 
Glegg’s x, y, z. Thus, in application of Liepmann’s spectrum, we set 

= Yo 

k. 2 - a (70) 

k 3 = v 

and recall that we have shifted the a integration in Equation 56 according to a=of+2nk/h. 

The normalized form for the Liepmann spectrum is 


(f>22 - 


w 2 21? L 2 (Kf + K^) 

Ustat nl + L2 ( K r + &2 + ^3 )] 3 


(71) 
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where I~A/R. In the Glegg routines, wavenumbers are normalized by chord c so that the 
parameter c = c / R appears in the formulas connecting the cascade wavenumbers and the 
turbulence wavenumbers: 


Aj — y 0 / c 

K 2 =atc =a' Ic+kBI Q 22 (72) 

A 3 = v/c 

where for A? , we have used 2 jvR/g - B and h/g~-Q 2 2 . Glegg’s wavenumbers are normalized 
by chord, i.e. y a = y Q c , etc. i; 

5.2 Axi-symmetric Spectrum 

When external turbulence is drawn into an engine inlet, it is distorted in the axial direction, 
producing elongated eddies at low forward speed and flattened eddies at high flight speed. For 
this kind of flow, Kerschen and Gliebe developed a model for axi-symmetric turbulence <ref n) . 
They gave the following general form for the axi-symmetric turbulence spectrum 

<Pij( k) = [k 2 Sjj - k t kj]F + [{k 2 - (k • A) 2 )8 ij - kjkj - k 2 AjAj + k • A(A t kj + k { Aj )JG (73) 

where A is a unit vector in the direction of symmetry and F and G are given below in terms of 
the wavenumber component k a aligned with the symmetry direction and the polar transverse 

wavenumber component k, = ^jkf v + kf z . 


and 


2 u 2 t / 4 

_ ~ u a l a l t 


7 

(i + el k 2 a + elkff 

Cl 

ii 

2 0 2 

2 U L--^L~\ 

«« >a J 

F 


(74) 


(75) 


The axial and transverse length scales i a and i t are independent as are the axial and transverse 

? 7 

intensities u a and uf . However, they are subject to the restriction 



(76) 


We assume the turbulence symmetry coordinates to have their x axis aligned with the flow in 
stationary coordinates and their z axis radial. To apply Equation 73 directly in Glegg’s cascade 
coordinates, we must establish the relation between symmetry coordinates of the turbulence and 
Glegg’s cascade coordinates. Transformation from one system to the other involves a sequence 
of rotations just like those shown in Figure 4 except that the starting coordinates are already 

rotated about the zj axis to the 6 S direction. Thus, the rotation matrix Q is the same as Q in 
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Equation 9 but with 9 replaced by AO - -( 9 S - 9) . Then, the transformation of wavenumber 
components is given by the inverse, i.e. 

k a = Q\\Y o + Q2\ a + Ql\ v 

k,y -Q\27 o + Ql2 a + 03 2 V (77) 

k tz - QnYo + Qn a + 033 v 

where 

0, , =cosAOcosy/ s + sinAtfsiny/'/ sin^ 9 
Q n = sinAflcos^. -cosA0siny// siny/ v 
013 = -COSyr/Sin^ 

021 = -sin A9cosy/[ 

022 = cosA^cosy / - 1 ('«) 

032 =-sin Wi 

0 3) = cosA#siny/ s . - sin A0siny^/ cosy/ 5 

032 = sinA^siny/^ + cos A0siny/7 cosy/ v 

033 = cos y/ j cosy/^. 

k-X in Equation 73 is most easily computed in the turbulence symmetry coordinates where 
X=(\, 0, 0 ). Then k-X. = k a . Otherwise in Equation 73, we need X in cascade coordinates, i. e. 


~x x r IT 1 ] 

X y = 0 x 0 = 0 21 (79) 

aJ L j L°J [a! 

Now, we apply Equation 73 directly in Glegg’s coordinates, where k t =y 0 , k 2 =a , and k 3 =v 

(f> 22 = (k 2 - a 2 )F + [(k 2 - kl) - a 2 - k 2 X 2 y + k a (2A y a)]G (80) 

where k 2 = y 2 0 + a 2 + v 2 . Substituting X y = 0 21 and re-arranging, we find 

<t>22 =(Xo + v 2 )F + K1-02 2 i)* 2 ~ k l -« 2 +2k a Q 2] a)]G (81) 

The non-dimensional forms for programming the axi-symmetric spectrum are 

<1>22 =(r 2 o + v 2 )F!c 2 + [(1 — 02 1 )7(. 2 — Kq - a 2 tc 2 +2K a Q 2] a/c)]G (82) 
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( 83 ) 


and 


with the requirement that 


_ u 2 i r 4 

L' _ J: “a _ _a 

* 2 (1 + iIkI + l]k}) 3 


G r = 


2 


ufj_U 

ulru 


2 

stat 

2 

stat 



IF 


2 f r t2 j2 

2 U LL— S J£L > 
w 2 / u 2 ~ I 2 

u a ' u slat lj a 

The normalized wavenumbers are 

K a = (0i ]X (; + 02i« + 03 i v )( ^ 
K/y = ( 012^0 + 022 ^ + 032*0 t C 
& tz = (013^0 + 023 a + 033^)/? 


and finally, 
and 


4 = 4 + 4 


K z = K 2 + Kf =(rl +a 2 + v 2 )/< 


(84) 


(85) 


( 86 ) 


(87) 

( 88 ) 
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SECTION 6 

APPLICATION OF THEORY & FORMULAS FOR CODE BBCASCADE 


To prepare the theoretical equations from Sections 4 and 5 for programming, we first write them 
in non-dimensional form and then modify them for spectrum level in decibels. Also, we address 
the radial integral over the turbulence field and limits for the wavenumber integral corresponding 
to cut on waves. 


6.1 Non-dimensional Forms 

To obtain power in any band of frequency, we should integrate the spectral density of 
Equation 64 over a range of a ) . For small enough bandwidth, however, we will simply multiply 
by A co and then multiply by 2 to account for negative frequencies. We call this integrated 
quantity PWR ± which, with the above definitions applied to Equation 64, becomes 

PWR ± = 2n2 — & Aa> f [ F ± (cc,v)Y \\<& 72 dz 0 ]tfti' (ft (89) 
j3s e (M + Q l2 M T )Q 23 JJ ^ L J 


where M^OR/a is tip (or local radius) rotational Mach number, and where the cascade power 
response function (in non-dimensional form and without the constants leading Equation 64) is 


F ± (a vt = - 

(i cD(4) 

2 

i yu , v ) 

l c 2 (K 2 e - fl) 


(90) 


We will do the radial integration over the turbulence (the inside integral) in a strip sense so that 
that integral can be approximated by 0 2 2 A R . T he result for one strip of width AR is 


2n BM slal 6) Ao) —\p g 2 RARa\ f f F ± (J «,P )£ 0 22 da’dv (91) 

Ps e (M + Q {2 M T )Q^° J JJ V 


where co = coR / a . The quantity in square brackets has the dimensions of power; the remainder 
of the expression is dimensionless. 

Power level in decibels is given by 


PWIr =10 log 10 


PWli 
PWR ref 


12 

where the power reference is 10' watts or 

PWRref = 0.73756 x 10 -,z ft lb / sec 
The result is that our working formula for the computer code is 


,-12 


(92) 


(93) 
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( 94 ) 


PW_ 

PWR rff 


7 7 __ 

2 JT BM; ta t co A a) 
fis e ( M + O] 2 M j ) (933 


\AP 0 RARa\ f f 

l'WR m , JJJ 


F ± (a,v)^ <t >22 da'dv 
k 


2 

We have used p n a =1.4/), where 1.4 is the ratio of specific heats for air and P„ is the 
ambient pressure. 

For use in the Glegg algorithms representing his D function, the non-dimensional 
chordwise wavenumber from Equation 65 


r, 


(<y- M z l’)-(Q22Gc' + Qt, 2 v )^T 


M + O Mn 


(95) 


..ere cb = coc/a is a non-dimensional frequency in the observer frame. The non-dimensional 
frequency (in cascade-fixed coordinates) from Equation 66 is 

— = o) — (Q\~>y a + Qtjo' + Q 22 F) Mt (96) 

a 

and for Glegg’ s shifted Helmholtz frequency to g 

® (co — W v)c A 

= — = (co - M z v) - ( Q\ 2 Y 0 + Qi2 a + Q 32 v )Mt (97) 

a a 

where M z ~W/a. Glegg’ s reduced frequency k (normalized by chord c) is 

^ 1/ c co K4 -rV A/ 7 

K = ~2 = P2 1 (012/o + 022 a + 032 v ) — 2 

a/ r p l p 2 


The remaining items to be dealt with for programming are the limits on the integrals and 
sum in Equation 94. These are treated in the following sections. 


6.2 Integration and Summation Limits 

The double integral in Equation 94 is over the a' , v wavenumber plane. There is an 
ellipse in that plane such that wavenumbers inside it correspond to cut on waves and 
wavenumbers outside correspond to cut off waves. In this section we find the equation for the 
ellipse and provide a scheme for determining loop limits (do-loops in the code) for the 2 
wavenumber integration variables. 

Cuton is governed by the square root in Equation 90. When the argument of the square root 

E = *l- fl (99) 

is positive, the waves are cut on and thus the cuton boundary is given by E = 0. We only need to 
treat k = 0, in which case Equations 5 and 6 from the Introduction give 
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( 100 ) 


K 


2 


_ (cr+KMd ) 1 _ n 

P 1 £ 


This is the ellipse in the cr, v wavenumber plane shown below in Figure 10. 



Figure 10. Sketch for discussion of wavenumber integration limits 

Solution for the a and v limits for the general rotor case involves algebra too tedious to 
present here. Instead, we present the derivation for the stator case, which is much simpler, and 
illustrates the concept. Then, the formulas for the general case are simply presented as 
documentation for BBCascade. 


Limits on Wavenumber Integral for Stator Analysis 

Equation 100 can be solved for cr as follows 

( 101 ) 
( 102 ) 


where, from Equations 3 and 5 


o = -KMd±s 0 ,\k l - v 


K - 


'P l 


co - vW 

ap 1 


We will use the outside loop for v and integrate over the range of a corresponding to <Ji < cr 
< &2 where 


cr, =-/cMd-s e J/c 2 - ^/ 2 


(103) 


-2 \P / 


(72 = -AT Md + s e 

as shown in the sketch. The extremes of v occur when the square roots above are equal to 0, 
i.e. for 

y. , = 0)1 a - (104) 

'• 2 W!a±P 
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or in non-dimensional form 


^ [}- M a sin y/ x 

p = 

2 /?+ AY 0 sin y/ s 


(105) 


v\here o.t = oxfa and M 0 = UJa . Finally, since <j~/od+ah the range of 
to Equation 1 03 is given by 





«2 = -(y 0 +ic M) f + 


/r 



« corresponding 


(106) 


This can be re-arranged and non-dimensionalized as follows: 


(107) 


a, =-(k/ M) tan z-^fi 2 * 2 -v 2 /cos z e 

a 2 = -(if/ M) tan z + yf S 2 k 2 - v 2 /cos 

where k -kc . These are the equations used for the stator analyses in Refs. 2 and 3 before the 
rotor theory had been completed. 

Limits on Wavenumber integral for General Case 

The limits for the rotor case were derived in separate notes and the results are presented here. 
Limits for vare given by 


v ],2 ~ 


CO 


( M z + M t Q 32 ) ± /f ^2MM r Q n - (Q} 2 + Ql 2 ) Mf 

and the a limits for any v in the above range are 

1 

an = — - 


i ~ ((oQ\2 + MQ 32 v) + 

M + Qj 2 M j* 

l 

^ ^2 A 2 - Mj J 

s e 


(- 2 - A } A 2 - coM T ) ± (coA] - M t A 2 ) 2 - P 2 v 2 A 2 + p 2 M 2 v 2 

s e V s e J 


(108) 


( 109 ) 


where 

and 


A, =P 2 -MM t Q ]2 

o _ G)„C 

^2 = P^Qn v 2 


( 110 ) 

(111) 
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These equations for the wavenumber limits are coded in BBCascade. They can be shown to 
reduce to the stator results above for M r = 0. Also, they have been tested numerically for 
correct behavior by one step in a or v and verifying that the argument of the square root in 
Equation 90 was negative. 

Summation Limits 

In the code the sum is initialized by the k = 0 term. Then, the k loop runs from 1 to a large 
number (30). The positive and negative k terms are added in pairs until a convergence criterion 
is satisfied. 


6.3 Sample Calculation - Comparison with ADP Test Data 

Figure 1 1 shows comparisons with scaled sound power spectral data from a model test in a 
wind tunnel at NASA-Lewis. Note that the noise below 300 Hertz has been traced to sources 
other than the fan. The stator was modeled with a single geometry and mean flow appropriate 
for the tip stator station. Turbulence scale and intensity were back-figured to provide a good 
match between the noise prediction and test data. The scale, A, was taken at 3.5% of the radius 
and the intensity was 2% of the mean flow velocity. Although the turbulence properties are 
represented simply by the Liepmann spectrum with a single intensity and length scale, the fit is 
satisfactory. (Issues of turbulence inhomogeneity were discussed above where it was shown that 
the theoretical equations for sound power presented herein apply equally well to the case of 
inhomogeneous flow, provided that the turbulence specrum is formally averaged in the gapwise 
direction.) The overall spectrum shape in Figure 1 1 is good and the high frequency slope 
matches that of the data. Furthermore, the split between the upstream and downstream power is 
excellent. 



Figure 1 1 . Comparison of theoretical and test noise spectra for scaled model data. 
Turbulence intensity and scale were chosen for a match to the noise data. 
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SECTION 7 


INTERPRETATION OF THE NOISE EQUATION - 
SPECTRA OF TURBULENCE AND CASCADE RESPONSE 


In the preceeding sections we derived equations giving sound power spectra (typified by l'igure 
1 1) in terms of the turbulence and cascade response spectra. In this section, before proceeding 
with the series of noise calculations in Section 8, we examine the behavior of the turbulence and 
cascade response spectra and how they work together to produce sound. To simplify this as 
much as possible, we consider the case for a stator with radial vanes. Under these conditions, 


Equation 56 becomes 


P s e k 




-/} 


dy 0 da d v 


( 112 ) 


where we have assumed that the turbulence is constant over a band ot radius A R and we are 
using V for vane count rather than B. The y 0 integral can be performed via the delta function 
with the result 


n ± Vp 0 n 2 c A co AR 

**0) ” o 11 
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JJ ^22 ( 


CO 

Tj, 


-,a, K 


>1 




da dv 


(113) 


The points below are best made in the context of circumferential duct modes. To simulate the 
duct, note that the kinematics of the excitation waves (at constant radius) would be expressed by 

e ‘( k xd x d +!> R yd) rat h er than the form actually used in the above analysis: e^ r ° x+ay ' > since the 
input field must be periodic on the interval 0 < y,j < 2n R . From the first exponential, we can 
see that the interblade phase angle (phase at y d = 2n R/V ) is 


<j = 


2nn 

IT 


(114) 


for the excitation waves. Then we can deduce, say from Equations B-6 and B-9 in Appendix B 
that the interblade phase angle of the response waves is 


2ir 

a-2nk = ~(n-kV) 
V 


(115) 


This is reminiscent of periodic rotor/stator interaction (Tyler-Sofrin) theory where wakes from 
the rotor with circumferential mode order nB scatter into circumferential order m=nB-kl at the 
stator. Since Glegg’s response potential depends on y 0 , a, and k only in the combination a- 
2xk,we can write 


a - 27ik = where m = n-kV (116) 

V 
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Then, in simplified form, the noise equation becomes 

00 00 00 

n ®= I I ( FH7t )« (H7) 

« ~ -00 A =-00 

where we have defined A' z - kT? ( the same definition as A^). Actual computed turbulence and 
acoustic response spectra at a mid frequency for Figure 1 1 (<y R/a = 50) are shown in Figure 12 
for k = 0. This shows that the sound spectrum is obtained as the “integral” in the transverse 
wavenumber ( n and K : ) plane of the product of the turbulence spectrum and the acoustic 
response spectrum. The turbulence spectrum is smooth and non-zero everywhere; the acoustic 
power response spectrum has a great deal of fine structure and is zero outside the cutoff ellipse. 
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Figure 12. Examples of the turbulence spectum and the acoustic power response spectrum (at it=0) 
for a mid-range sound frequency (w R/a = 50) 
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Figure 13 shows turbulence spectra at the top and acoustic response spectra at the bottom 
for higher and lower frequencies. It can be seen that both spectra broaden in the wavenumbei 
plane with increasing frequency and that the acoustic response spectrum maintains its fine 
structure independent of frequency. 



Figure 13. Turbulence spectra (top row) and acoustic spectra (bottom row) at 3 sound frequencies: 
coR/a = 10 (left) , 50 (center) . and 100 (riqht). 
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Figure 14 at the bottom shows the cascade respnse spectrum with several values of scattering 
index included. It can be seen that the acoustic response for each k is identical but shifted by 
the vane count V ( 30 in this example) through the mode formula m^n-kV At low frequency 
= 10, the sections of the response spectrum do not overlap. Comparison with the turbulence 
spectrum above in the figure shows that only small portions of the wavenumber plane contribute 
to the sound. Specifically, response to radial wavenumbers is limited by cutoff approximately to 
the range ±12 where the turbulence spectrum is near its peak. However, the response to 
circumferential wavenumbers is not limited by the acoustic response spectrum but by the rolloff 
of the turbulence spectrum. 

The right hand part of Figure 14 applies to the mid-frequency point of Figure 11. In 
accordance w ; *h Equation 100, the “radius” of the cuton ellipse increases but the sections of the 
acoustic response sinti only by V (independent of frequency), so that they overlap. Thus, some 
turbulence wavenumbers produce noise in multiple circumferential modes. 




Mid Frequency: 03= 50 







SBl 

■ 

II 

to 

u> 


o 50 100 


Acoustic 
Power Response 
Spectra 


I (PWRC.-tv 

Jt = -oo 


CIRCUMFERENTIAL WAVENUMBER n 


Figure 14. Contour plots of the turbulence and acoustic spectra for a several values of the scattering 
index k. 
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SECTION 8 

COMPARISONS WITH TEST DATA & RELATED PARAMETRIC STUDIES 


In this section we explore the capabilities of BBCascade in prediction of broadband noise of 
stators and rotors and establish sensitivity to design and flow parameters. Fans are represented 
via a single representative radius calculation rather than making any attempt to integrate over fan 
radius as is done in the BFaNS code. Thus, we do not expect perfect agreement but look for 
correct spectrum shapes and correct trends and approximate absolute levels. Three data sets are 
used: 

ADP Fan 1 (Section 8.1) 

22” low speed, adjustable pitch fan 
Tested in the NASA-Glenn 9x1 5 tunnel 

Allison Fan (Section 8.2) 

22” fan with combinations of stator lean and sweep 
Tested in the NASA-Glenn 9x15 tunnel 

Boeing Fan (Section 8.3) 

1 8” fan with variable geometry for noise source separation 
Tested in the Boeing LSAF (Low Speed Aeroacoustic Facility) 

In each of the subsections below (Sections 8.1, 8.2, 8.3), the format is to check predictions for a 
few data points to verify correct behavior of the theory and to establish base cases for further 
studies. Then we present parametric variations on the base cases to provide an understanding of 
the noise generation process. 


8.1 Variations on ADP Fan 1 - Turbulence Scale, Mach Number, 
Stagger and Chord 

Since noise predicted by Equation 94 depends on only a small number of input variables, it was 
straightforward to compute the variation of noise with turbulence properties (intensity and scale), 
mean flow (Mach number) , and geometry (stagger, vane count, and gap). The base case for the 
variations was essentially the same as that for Figure 1 1 except that the vane count was changed 
to 30. Stagger was 30°, gap/chord was 0.8, Mach number was 0.5, turbulence level (upwash 
component) was 2% of the mean flow velocity and the integral scale was 3.5% of the fan radius. 

First, recall how the curve fit in Figure 1 1 was obtained. Mach number, stagger, chord, 
solidity, radius, and area were fixed leaving only turbulence instensity and scale as free 
parameters. From Equation 94 it is evident that spectra simply move up and down with 
turbulence intensity (3 dB per doubling of turbulence energy). Turbulence scale, however, is 
more interesting, as shown in Figure 15. (In this and succeeding plots, the solid curve is the 
base, or reference case.) At high frequencies, the level varies as 1/A, which could be deduced 
from the form of the turbulence spectrum in Equation 69. At high frequency, the 1 in the 
denominator becomes insignificant and all but one of the A’s cancel. At low frequencies, the 
effect reverses and scale has a powerful effect. Thus, to achieve the match in Figure 1 1, intensity 
was varied to move the theory curves up and down and scale was varied to change the balance 
between low and high frequencies. 
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Figure 16 shows the variation of noise with Mach number at constant percent turbulence, i.e. 

o I o 

for constant w /U stat in Equation 71. On the upstream side, the curves generally move up in 
level and frequency with increasing Mach number. Behavior on the downstream side is similar 
except that the spectrum tends to flatten at the higher Mach numbers. Peak levels of these curves 
were cross plotted against Mach number in Figure 17 to find if a simple power law ( A/" ) would 
describe the behavior. All curves were referenced to the M=0.5 values. On the downstream side, 
M 5 fits reasonably well. On the upstream side, no exponent works perfectly but 4.5 fits the 
points between M— 0.5 and 0.7. 

Figure 18 shows variations with cascade stagger angle. The effect is weak on the 
downstream side and almost zero on the upstream side. 

Finally, we look at variations with vane count and gap/chord ratio. Vane count is a major 
consideration for harmonic noise; large numbers lead to cutoff (exponential decay in the axial 
direction) of the lower harmonics. In designing stators, there is usually a preferred solidity 
(inverse gap/chord ratio) for performance purposes but some latitude for vane count at constant 
solidity. Figure 19 shows broadband noise variation with vane count at constant solidity. At the 
mid and high frequencies the noise is roughly proportional to vane number (PWL«10 log V). At 
lower frequencies there is little effect. Thus, we are driven in different directions for tone and 
broadband noise. Figure 20 shows that variations in gap/chord ratio have almost no effect for 
constant vane count. This was a surprise since acoustic response for individual input waves is 
very sensitive to gap/chord; however, when the turbulence is integrated over a broad range of 
wavenumbers, the sensitivity to gap/chord is much weaker. 
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Figure 15. Variation of stator noise with turbulence scale. L/R is turbulence integral scale/fan radius. 
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Figure 16. Variation of stator noise with inflow Mach number at constant percent turbulence. 
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Figure 17. Variation of spectrum peak levels with Mach number. 


Curve fits represent PWL-A/ 1 
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Figure 19. Variation of stator noise with vane count at constant gap/chord ratio 
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8.2 Allison Fan and the Effects of Lean and Sweep on Stator Noise 

For another check of the theory, we compared calculations with the lean and sweep data 
from the Allison fan ref 12 shown in Figure 21 . The stator at the top in the figure has its tips swept 
axially; i.e. per the nomenclature of this report, it is swept 30° in the duct system ( y/ v — 30 ). 
At the bottom Figure 21 shows another stator with both lean and sweep in the duct system 

= 30° and y/ t = 30°). Allison and NASA tested both of these stators plus radial stators at 
2 axial positions. One axial position was at the root position of the swept stator and the other 
was at the tip position of the swept stator. 

Figure 22 shows sound pressure spectra in the aft quadrant for the 4 stator configurations. 
Although the figure is difficult to decipher, it can be seen that sweep eliminates the tones. 
Furthermore, sweep provides broadband noise reduction over a wide range of frequency and lean 
added to sweep provides a further improvement. Reference 12 gives the vane count and the 
sweep and lean angles; however educated guesses had to made for the other parameters needed 
for spectrum calculations. We used the same non-dimensional parameters as for Figure 1 1 
except that the Mach number was reduced to 0.3 (based on Allison’s 50% design tip speed). 
Figure 22 at the bottom shows that the theory follows the Allison data very well. The spectrum 
shapes are good and the reductions for sweep and sweep+lean match the data as well as can be 
determined from the curves. 

The remainder of this section explores theoretical trends with lean and sweep and explains 
the Allison results. 

Lean and Sweep Parametric Studies 

In this section we examine trends in noise reduction with various combinations ol lean and 
sweep. This is done in the context of the 2 conventions for defining lean and sweep: the 
“cascade system” of Figure 4 and the “duct system” of Figure 5. 

The following spectra were all computed using the same geometry, mean flow, and 
turbulence properties used for Figure 11. Only the downstream sound power spectra are shown 
since they are higher and look similar to the upstream spectra. Figure 23 shows the effect of lean 
and sweep for angles defined in the duct system. Lean has almost no effect except to move the 
dip in the spectrum. (In reference 3 it is shown that the frequency at this dip corresponds to a 
half wavelength fitting between 2 vanes. Leaning the vanes changes the normal gap - Glegg’s 
h.) Sweep is more effective; the slight increase at low frequency would not be seen in fan test 
data because of other sources. Figure 24 shows lean and sweep effects for angles defined in the 
cascade system. Here, again, lean is very weak but the sweep effect is stronger than in the duct 

system. .... . , 

The previous curves were for lean alone or sweep alone. Since it is important to know how 

the 2 effects work together, combinations of lean and sweep were run and the peak values from 
the downstream spectra were plotted in Figure 25. The top of Figure 25 applies to the cascade 
system (sweep in the plane of the vane). It shows noise reduction versus sweep for lean angles 
of 0° , 1 5°, 30°, and 45°. The curves show that sweep is the dominant effect. In this system 
sweep reduces the component of Mach number normal to the leading edge; lean does not. The 
solid curve shown in the figure 


AdB = -30 log] 0 (cos y^) 

provides a reasonable fit to the theory. (Recall that the lengthening of the leading edge is 
included in the theory. Without this, the factor would be 40 rather than 30.) Negative lean 
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angles are not shown. However, the curves apply to negative lean if the signs of the sweep 
angles are reversed. 

The bottom of Figure 25 applies to the duct system, i.e. lean and sweep angles defined via 
the front and side views of the fan. Here the picture is complicated due to a strong interaction 

between lean and sweep. However, -30 logio(cos^) is still valid (no hat on y/ s ). We need only 
apply Equations 15, which relate angles in the duct and cascade systems, to convert the solid 
curve from the top of Figure 25 to the multiple curves at the bottom. The curves track the 
theoretical peaks reasonably well and can be used for design guidance. 

, mally. refer again to the Allison data in Figure 22. Sweep and lean were both 30°. 
Calculated results shown in Figure 25 predict that 30° sweep by itself produces about 1.1 dB 
benefit and that 30° lean added to the sweep brings the total reduction to 3.9dB. This is in 
approximate agreement with the Allison data. If they had leaned in the opposite direction, the 
i. i eTect vv^-ild hav'e been near 0 dB. The curve-fit plots in Figure 25(the solid curves) indicate 
slightly different results: 1.5 dB for sweep and a total of 3.0 dB. Note that is was not luck that 
caused Allison to lean in the more favorable direction; this is the direction preferred for tone 
noise reduction. 

The curve fit without the BBCascade calculations is shown in Figure 26 with sketches 
defining the lean and sweep conventions. Note that this is simply a plot of -30 logio(cos^) in 
the duct (or “hat” ) system for 30° stagger. Details will be different for different stagger. 
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C-96-4831 


Photograph of the partially-assembled fan 
stage showing the swept-only stator. 



Downstream view of the swept + leaned stator. 

Figure 21. Allison fan disassembled to show swept stator (top) and swept and leaned stator (bottom). 
Reproduced from Woodward, ef a/, ref. 12. 
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Figure 22. Top -variation in sound pressure spectra with lean, sweep, and axial position of stator. 
Allison fan at 50% speed/1 26° from inlet axis. Bottom - predictions for aft sound power spectra. 
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Figure 23. Variation of stator noise with lean and sweep angles as defined in duct coordinates 
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Figure 25. Variation of stator noise with combinations of lean and sweep. Top: angles defined in 
cascade system. Bottom: angles defined in duct system. Solid curves are AdB=-30log 10 (cosv s ). 
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Figure 26. Plot of noise peak noise reduction from curve fit: AdB=-30log 10 (cosvj/ s ) . (Same as in 
proceeding figure. 
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8.3 Boeing Fan Analysis and Variations 

This section checks BBCascade against Boeing’s data for inlet turbulence/rotor noise, rotor wake 
/stator interaction noise, and mode order/frequency plots. 

Inlet Turbulence / Rotor Noise 

The Boeing 18” noise research fan was designed for source separation by experiment. It can be 
run without stators and with the duct boundary layer upstream of the rotor removed. By 
comparing rotor noise with and without the boundary layer, the contribution of boundary layer 
turbulence/rotor noise can be determined. Then by comparing noise without the boundary layer 
with and without stator, the contribution of the stator can be determined. The results of this 
source separation are shown in Figure 27 (reproduced from ref. 1) for inlet sound power spectra 
at 70% rpm and aft sound power at 55% rpm. (100% rpm was equivalent to 1300 ft/sec 
rotational tip speed.) The boundary layer source is not dominant but it can be seen to affect the 
total level in both cases over a broad range of frequency. The bold curves in Figure 28 are a test 
of BBCascade’ s rotor noise capability. These were computed using the geometry and mean flow 
conditions at the rotor tip. For turbulence intensities, we used the values at a distance 50% of the 
boundary layer thickness (8) from the wall were used (4% / 3% for the axial / transverse 
components). For turbulence axial and transverse scales, we simply adjusted them to produce a 
good fit to the noise spectrum as shown at the bottom in Figure 28. The values that worked best 
were axial / transverse scale = 7% / 2% of fan radius. This is not very satisfactory since the axial 
scale of the streamwise component that would be deduced from ref. 1 at 50% of 8 is roughly 
3.5% of fan radius. We do not have an explanation for this discrepancy. Since the transverse 
scale was not measured, we simply comment that 2% of fan radius seems reasonable. 

For the aft noise, at the bottom of Figure 28, the fit is very satisfactory. The spectrum peaks 
are centered at multiples of blade passing frequency. Agreement for the inlet noise is not as 
good; however, the Boeing report points out that there is a rotating instability that influences the 
spectrum. Of course, this is not accounted for in the theory. 

The theoretical curve from the top of Figure 28 was used as a base case to study effects of 
changing turbulence length scales. Axial scale was 7% of fan tip radius and transverse scale was 
2% of tip radius. Figure 29 shows variation with transverse scale in 2:1 steps. Note that this 
parameter is a very strong driver and must be known with some accuracy for acceptable noise 
estimates. Figure 30 shows variations with axial turbulence scale again in 2:1 steps. Axial scale 
does not change overall level strongly but instead controls the localization of energy around the 
blade passing frequency peaks. It can be seen that the 3.5% axial scale that would be deduced 
from the Boeing hot wire data produces a spectrum with peak-to-valley ratios that are too small. 

Usually, we think of fan noise spectra containing pure harmonics (peaks with zero 
bandwidth at multiples of blade passing frequency) and a smooth broadband component. 
However, the results just shown are in between harmonic and broadband. The peaks are 
centered at the BPF harmonics but have finite bandwidth. This is caused by the rotor interfering 
with turbulent eddies that are long enough to chopped by more than one rotor blade. A pure tone 
would be produced by an infinite number of chops but narrow band peaks are produced by finite 
numbers of chops. The first good model for this effect was due to Mani (ref. 13) based on 
isotropic turbulence. The present author brought in the essential effect of anisotropy in 1974 
(ref. 14) based on a cruder model of turbulence interaction. It is interesting to see how long the 
axial scale must be to produce the spectrum peakiness. The sketch on the next page shows an 
eddy interfering with a rotor with blade gap g and stagger 6 . The axial distance between 
chops is X — g cot 0 and the gap is g-InRJB. Combining these formulas with the axial scale- 
on-radius used for Figure 28, A/R = 0.07 and a stagger of 60° leads to the ratio of axial scale to 
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axial distance A/A, — 0.4 . This is surprising since it indicates spectrum peaking occurs when the 
axial scale is shorter than the axial gap. Since Boeing obtained both noise and turbulence 
spectrum data supporting this result, we must conclude that it is not the mean eddies that are 
dominating the noise but the eddies that are larger than average. 

Figure 31 shows the effect of changing blade count on turbulence/rotor interaction noise. It 
can be seen that doubling blade number cuts the number of peaks in half and raises their level. 
Also, the peak to valley ratio increases. Finally, 

Figure 32 shows the effect of sweeping the rotor 
tip 45° (in cascade system, i.e. sweep in the plane 
of the rotor blade). There is a significant benefit 
but only at the higher frequencies. 

Rotor Wake I Stator Interaction Noise 

For the Boeing tan aata, we also look briefly at 
trends in stator noise with rotor tip speed and with 
vane count. Figure 33, copied from the Boeing 
report, shows inlet and aft spectra for 3 rotor 
speeds and 3 vane counts (15, 30, & 60). Noise 
merci es with tip speed, as expected. Noise also 
increases with vane count; in fact, for the lower 
speeds, sound power is proportional to vane count 
(3 dB per doubling of vane count). At the higher 
speeds, the inlet noise results are more difficult to 
characterize, perhaps due to rotor transmission 
effects. Note that the solidity for the 1 5 and 30 
vane stator were the same but that the solidity for 
the 60 vane stator doubled. 

The corresponding theoretical curves are 
shown in Figure 34. To produce this figure, the geometry and mean flow from the 85% radius 
were matched for the 3 rpms and then turbulence intensity and scale were adjusted to fit the 
middle noise spectrum in the aft quadrant (78% rpm, 30 vanes). The other vane counts were 
computed assuming that the flow properties were unchanged. Considering that the fan is 
represented by a single representative radius, the trends with rpm and vane count are predicted 
very well. In particular, over most of the spectrum, broadband noise is proportional to vane 
count, an important result for design considerations. 

Mode Order Frequency Plots 

Finally for the Boeing data, we examine mode order/frequency 
plots to be sure that BBCascade has the correct behavior at that 
level of detail. Boeing used rings of Kulite pressure transducers 
in their 18 inch research fan to determine the distribution over 
circumferential mode order (of pressure at the outer wall) as a 
function of frequency. Examples are shown in Figures 35 and 37, 
which were copied from the Boeing report (ref. 1). Figure 35 is 
from a ring of Kulites between rotor and stator and is believed to 
be dominated by stator noise. The V-shaped pattern is the cuton 
limit, although there are pressure modes outside this limit because 
of the near field location of the transducers. Co-rotating modes 
(m > 0 ) are dominant upstream of the stator. This is expected, based on the sketch at the right, 
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since waves propagating normal to the airfoil surfaces should couple most efficiently. The 
“cuton V” is skewed because the swirl between rotor and stator cuts on counter rotating modes 
at lower frequency than co-rotating modes. Figure 36 shows predicted mode order/frequency 
plots for sound power both up and downstream of a stator. On the upstream side (top portion of 
figure), the co-rotating modes clearly dominate and the cuton V is correctly skewed since the 
case was run with swirling input. We don’t expect exact agreement with data since the 
experimental plots are for pressure modes and the theory is for sound power; however the correct 
qualitative behavior has been verified for the stator. Figure 36 also predicts downstream of the 
stator that the counter-rotating modes will dominate, although Boeing made no measurements at 
this location. Note that the cuton V is not skewed because the case was run without swirl. 

What behavior do we expect upstream of the rotor? Based on a sketch analogous to the one 
above, we might expect that counter-rotating modes would dominate since rotor stagger has the 
opposite sign. However, the experimental results from Boeing indicate that again the dominant 
modes are co-rotating. (This issue is clouded somewhat since the data include noise from 
trailing edge noise as well as duct turbulence interaction.) However, the theory, shown in 
Figure 38 predicts dominant energy in the co-rotating modes also. Apparently, the motion of the 
rotor biases the energy distribution in this direction. The theory must account for this correctly 
via the frequency scattering, which is different for each mode order. 
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Figure 27. Boeing’s source separation results (ref. 1) showing experimentally derived rotor alone 
noise (no stator, no duct boundary layer), duct turbulence/rotor interaction noise, and rotor wake/ 
stator interaction noise. 
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Figure 28. Calculated results from BBCascade for duct turbulence/rotor interaction superposed on 
Boeing’s results from preceeding figure. 
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Figure 29. Variation of rotor noise with changes in transverse scale of inflow turbulence. Middle case 
is the match to Boeing data shown in Figure 28. 



Figure 30. Variation of rotor noise with changes in axial scale of inflow turbulence. Middle case is the 
match to Boeing data shown in Figure 28. 
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Figure 31 . Variation in rotor noise with changes in blade count (at constant gap/chord). 



Figure 32. Variation of rotor noise with sweep. 45° sweep as defined in cascade system. 
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Figure 33. Figure from Boeing report (ref. 1) showing net stator noise at several fan speeds and with 
vane counts of 15, 30, and 60. 
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Figure 34. Calculations from BBCascade approximating conditions for Boeing fan in 
preceeding figure. 
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Figure 35. Mode order/frequency plot from Boeing report (ref. 1) from Kulite array between rotor and 
stator. Note dominance of co-rotating modes (m > 0). Also, note that contra-rotating modes ( m < 0) 
cut on at lower frequency than co-rotating modes. 
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Figure 36. Mode order/frequency plots from BBCascade for stator showing distribution of sound 
power over circumferential order (summed over all radial wavenumbers). 
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Figure 37. Mode order/frequency plot from Boeing 
rotor. Test configuration without stator. 


report (ref. 1 ) for transducer array upstream 


of 
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Figure 38. Mode order/frequency plots from BBCascade for rotor showing distribution of sound power 
over circumferential order (summed over all radial wavenumbers). 
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SECTION 9 

CONCLUDING REMARKS 


A theoretical formulation was developed for broadband noise generated by turbulence 
convected into a cascade. The theory was derived by applying Glegg’s harmonic recitlincar 
cascade theory to the problem of random inflow with resulting equations giving upstream and 
dowstream sound power spectra as integrals over the 3D turbulence spectrum. Within the 
limitations of the rectilinear cascade model, the analysis was made as general as possible. In 
particular, it treats noise from rotors as well as stators; the blades or vanes can be leaned and/or 
swept; and turbulence that can be inhomogeneous and can have a limited form of anisotropy. 
Two spectra for turbulence are provided: one is the Liepmann spectrum for homogeneous, 
isotropic flow; the other for axisymmetric turbulence with stretching in the direction of the mean 
flow. Operation of the theory was explained in some detail with actual plots of turbulence 
spectra and cascade power response spectra. 

The theory is embodied in code BBCascade. This was developed for application to turbotan 
noise but does not truly provide a fan prediction method since the mean flow and geometry are 
independent of spanwise position. For the fan application, BBCascade has been incorporated in 
code BFaNS (described in a reference) with integration over varying geometry and flow 
properties in the spanwise direction. BFaNS also accounts for trailing edge sources at the rotor 
and stator. 

With the constant geometry and mean flow it was instructive to perform a series of 
parametric studies to explore the fundamental behavior of predicted noise. The format for these 
studies was to first compare predictions with data from a model fan, then find flow parameters 
that give a convincing prediction of the noise, and finally perform parametric variation on the 
geometry, mean flow, and turbulence. This provides a basic understanding of the methodology 
and of the sensitivity of noise to various design and flow parameters. Fan models included a 
Pratt & Whitney 22 inch model fan, an Allison 22 inch model fan with swept and leaned stators, 
and a Boeing 1 8 inch fan designed to separate the various broadband noise sources. 

Agreement with data is very convincing and explains several observed phenomena. Of most 
interest today perhaps, the reduction of broadband noise by leaning and sweeping stators is 
thoroughly investigated and the critical parameters are identified. Predicted noise is reduced by 
sweeping the vane in its own plane according to AdB = -301og 10 (cos^/ 5 ) , where y/ s is the 
sweep angle. Rotor noise also can benefit from sweep. 

Other items studied include sensitivity to geometry variations such as solidity, vane count, 
and stagger as well as sensitivity to mean flow Mach number. Predicted noise is roughly 
proportional to vane count in the important frequency ranges, i.e. PWL <x 101og 10 V where V is 
vane count. Accurate prediction of broadband noise obviously requires accurate specification of 
the turbulence properties. Noise increases uniformly at all frequencies with turbulence intensity, 

i.e. 3 dB per doubling of turbulence intensity, w 2 /u 2 . However, it is also shown that 
knowledge of turbulence scale can be even a more important driver: doubling or halving 
transverse scale can cause noise changes greater that 3 dB. 
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APPENDIX A 

CASCADE GEOMETRY FOR LEAN AND SWEEP ANALYSIS 


This appendix derives geometrical relations between duct coordinates (where chord, gap, and 
stagger angle are Cd, g , #) and Glegg's cascade coordinates coordinates (where chord, gap, and 
stagger angle are c, s, % ). In the process, we develop the coordinate translormation equations 
and also define lean and sweep according to 2 different conventions. Geometry for the 
derivation is presented in Figure 4 of the main text. 

Stagger, lean, and sweep are defined by 3 successive rotations about coordinate axes as described 
in Section 3. The rotations are non-commutative so that the order of the rotations must be part of 
the convention. 


Coordinate Rotations 

Stagger is the rotation from the duct 
system (xd, yd, z d) to the (x 0l y 0 , z 0 ) 
coordinates. The latter system is used in 
the main text for analysis of the 
turbulence Fourier transforms. The 
sketch at the right represents a 
cylindrical cut of the stator unwrapped 
onto a plane. Gap and stagger angle are 
g and 9. 



The transformation can be written in matrix form as 


*o" 

To 

II 

JO 

X 

i 

"CJ JJ* 

_ z o_ 


l z d\ 


where 



cos 9 
-sin# 
0 


sin# 

cos# 

0 


0 

0 

1 


(A-l) 


(A-2) 
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Lean is a rotation through angle y\ about the common x Q , xj 
axes into the primed system as shown at the right. Consider 
the vanes to have “hooks” on their leading edges attaching 
them to the yj axis. Thus, the cascade does not rotate as a 
unit. Rather, each vane rotates on its hook, maintaining the 
original gap g at constant radius. However, the normal 
gap, used in the cascade theory varies with lean per formulas 
given later in this appendix. 

The transformation is 



I 

* X 

1 

X 

§ 

11 

J 

G O 

1 

1 

1 


1 

G 

1 


(A-3) 


where 



0 

COS l// 1 

sin {// 1 


0 

-sin y / 1 

COS !)/ 1 


(A-4) 


Sweep is a rotation of the vane (in its own plane ) about the y ' 
axis through angle ig s ■ This establishes the chord at constant 
radius in duct coordinates as cj and the chord, measured 
normal to the leading edge as c with the relation c - Cd cos if/ s . 


The transformation is 



X 


x’ 

y 

X 

O? 

II 

y 

z 


z' 


where 


(A-5) 



cosy^ 

0 


sin 


0 -sin^ 5 

1 0 

0 cos y/s 


(A-6) 
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The net result of the 3 rotations is 


[sHaMaM&] 


(A-7) 


or 



cos # cos ^ + sin 6* sin y/ / sin 

-sm9cosy/ / 

cos 9 sin y/ s - sin 0 sin y// cosy/,. 


sin 9 cosy/ s - costfsin y / 1 sin y/ s 
cos 9 cos y / 1 

sin 9 sin y/ s + cos6>siny//Cosy/ 5 


-cos y/[ siny/ ? 

-s'my/j 
cos y / 1 cos y/ s 

(A-8) 


giving the desired transformation of the duct coordinates into the cascade coordinates. 


X 


x d 

y 

=[£]« 

yd 

z 


z d _ 


(A-9) 


The inverse transformation is also useful and is just the transpose of \Q\. Define the inverse as 


so that 


bHer' 


(A- 10) 


r i 
Xd 


X 

yd 

= [q]x 

y 

z d _ 


z 


(A-l 1) 


where 



cos 9 cos y/ s + sin 9 sin y/ / sin y/ s 
sXndcosy/ s -cos9s\ny / 1 siny/ 5 
-cos y/{ siny/ 5 


-smOcasy/i cost?siny/ 5 -sin#siny// cosy/ s 
cos9cosy/i sin#siny/ s + cos#siny// cosy/ s 
-sin y/f cosy// cosy/ 5 

(A- 12) 


Lean and Sweep as Viewed Parallel and Perpendicular to Fan Axis 

Another scheme to define lean and sweep is according to the front and side views of the 
unwrapped blade row, as sketched in Figure 5. From the front (looking along the Xd axis) the 
lean angle is y/ / and looking from the side the sweep angle is y/ s . In cascade coordinates, the 
leading edge of the reference blade is located at x=0, y=Q. Equation A-l 1 transforms this to the 
duct system: 

Xd=? 13 Z 

yd =c l23 z (A- 13) 

z d = ^33 z 
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From the front, the slope is 


and from the side, the slope is 


tan y/ s 




or 



^ - _ (J32 _ sin Os\ny/ s + cos 9 sin y/j cosy/ s 

033 cos y/, cos if/ s 

* _ 03i _ cost?sin^ y — sin^sin^/cos^ 5 
tan (// ^ — — * — - — 1 

033 cos y / 1 cos y/ s 


These equations can be inverted to give 


tany/7 = cost? tan y/f -sin^tany^ 

tan y/ s = cos y/[ (sin 6 tan y/i +cos6 , tani^ v ) 


(A-14) 


(A- 15) 


which must be solved sequentially. 

Other expressions required in the text can be obtained from sin 2 y/ s + cos 2 y/ s = 1 : 

sin^ = T A'", ' cos ^= / -f 33 0 (A- 16) 

V 031 + 033 3/031 +033 

Finally, for Equation 51, we need an expression for g ! {scosy/ s ) . For this, first note that the 
square root in the expression for cosy/ s above can be written 

V 03 1 + 033 = a/031913 + 033933 (A- 17) 

where q is the inverse (and therefore transponse) of Q . And since Qxq = 1, it can be seen that 
031^13 + 032923 + 033933 = 1- Therefore 

V 031 + 033 = V 1- 032923 
Similarly, from Equation A-23 (below) 

~ = V 0T 2 + 022 = a/921012 +922022 = a/ 1 “923032 
& 

Then, combining, Equations 16, 18, and 19, we find the desired result for Equation 51 

g _ i 1 

s cos y/ s 033 cos y / 1 cos y/ s 


(A- 18) 
(A- 19) 

(A-20) 
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Relations Between Definitions of Chord, Gap, and Stagger 

Given the geometry in the duct system in terms of chord c d , gap g , and stagger angle 0 , we 
can now find the corresponding quantities in the cascade (Glegg’s) system. To do this, we note 
that the coordinates of one point on vane #1 in duct coordinates are x d = 0, y d = g, z d = 0. From 
Equations B-l 1 and B-12, the corresponding point on vane #1 in the cascade coordinates is at 

x = Q \ 2 g = (sin 0cosy / s - cosflsin y / / sin y/ s )g 
y = Qn g = (cos0cosy / i)g (A-21) 

z - Q 32 g = (sin 0sin y/ s + cosflsin y/ t co sy/ s )g 


Since the leading edge is normal to the x, y plane in cascade coordinates, x and y are Glegg s d 
and h in the sketch below. 



c 


Thus, we have 

d = Q\2 g = {sm 0 cosy/ s -cos 9 s\ny/i siny/ 5 )g 
h = Ql2g = {cosO cosy/ i)g 


And, from these, we immediately write down the expressions for Glegg’s gap and stagger angle: 


s = ^d 2 +h 2 = gyjQ\2 + Q22 

(A-23) 

tan 

(A-24) 

We have already noted that Glegg’s chord is 


c = c d cos^ 5 

(A-25) 


Also, note from the sketch above that the velocity components are 

U — U 0 cosy/ s 
W = U 0 s\ny/ s 


(A-26) 


where U 0 is the full velocity into the cascade (assumed to have no radial component in the duct 
system). 
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Thus, for a typical application, we might be given chord c d , gap g, stagger 9 , lean and sweep 
angles, y/ / and y/ v in the duct system. Lean and sweep would be converted to the “no hat” 

system via Equations A- 15. Then the chord c , gap 5, and stagger x f° r Glegg’s geometry 
would be obtained from Equations A-25, A-23, and A-24 and his velocity components from 
Equations A-26. Interblade phase angle and spanwise wavenumber are discussed in the main 
text. 


Rotation from the “Duct” System to the “Hat” Coordinate System 

One more coordinate system is needed in the main body of the report for integration of intensity 
over the cascade face to obtain sound power. This is the “hat” system reached by rotating about 
the transverse duct axis yj through the angle \j/ s per the sketch below. The transformation 

from the duct system to the “hat” system is 

z d W s 

x d cos^ ? 0 sin^j x 

y d = o 1 0 X y (A-27) 

z d -sinyLj 0 cos^ 5 £ 

and hence the net transformation from the “hat” system to 
$ GI egg’s cascade system, namely 



is given by the matrix 

Q 





COS0 5 

0 

smy/ s 

X 

0 

1 

0 


-siny/j 

0 

cos^ 5 


or 

(01 i cos ^ - 0j 3 sin ^ s ) 0 12 (0n sin^ + 0 13 cosy^) 

Q= (Q 2 \^°Ws-Q 23 s[n ¥ s ) Q22 (02isin^ + 02 3 cosy>. ? ) 

_(03 1 cos ¥ s~ 033 sin (j/ s ) 032 (031 sin^ s + 0j 3 cos^j) 

and, by use of Equations A- 16 and Qxq = 1, an alternate form can be found 


(A-28) 


(A-29) 


(A-30) 
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(A-31) 


0 = 


0 n 03 3-01303 i 

yjofi + Qh 

02 i 03 .3- 02 3 03 i 
+ 033 


0 


02 

022 

032 


-032012 

^ 3 \ + 033 
-032022 

V 031 +033 
1 ~ 032 

^^ 0 ^ + 033 


This is used in Section 4 of the main text and in Appendix B. 
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APPENDIX B 

WAVENUMBERS IN VARIOUS COORDINATE SYSTEMS 


In the main body of this report, wavenumber vectors must be expressed in a variety of coordinate 
systems. This appendix addresses that issue in general and also derives a relation involving 
wavenumbers that is needed for Equation 42 in derivation of the acoustic equations. 


The first issue is easily addressed. In Appendix A several coordinate system were defined by 
rotations about coordinate axes. Matrices representing these rotations were derived and given 
names of Q with various subscripts, tildes, and hats. A basic theorem of tensor analysis is that 
the matrix that rotates position vectors also rotates other vectors. Hence, if 


then 


X 

y 



Q 

X 

X 

y 

z 


- 


z 

i 


Q 

X 

Kd 

kyd 

V 




Kd 


(B-l) 


(B-2) 


The other topic of this appendix is to re-express the exponential in Equation 42 in a form more 
suitable for integration. We extract the term 

= -4* -K4 7 + + •* + <»*' < B - 3 > 

n h 

from Equation 42. This expresses the kinematics of the cascade response waves. We will 
express this in the x,y,z system, identify the wavenumbers of the excitation waves in this 
system, and then re-write the exponential in the desired form. If we use Equation B-l to 
transform the x,y,z coordinates, E R becomes 

4" A ~e d fT 2 nk Cl Cl i a 

E r = [-4e n +(AT- + — )02l + ^31^ 

+[-*kQl2 +('L ± k~ + j- )Q22 + v Ql2]y (B-4) 

+[-4a 3 +<4 T + ^T^>2!3 + ^ 

n n 

The terms involving the T’s on the 2 nd and 3 rd lines cancel, leaving 

„ , „ + d <r-27lk.j\ A - 

E r = [-h k Q\\ +(A k — + - )Q2\ + v Qi\\ x 

h h (B-5) 

<j — 2/rk _ ^ g — 2nk ^ * a ~ 

+[( ) + K?32 + [( 7 )£?23 + V '0J3 ] z 

g h 
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where we have used Qi2~Qi2> Q 22 = Q 22 ■> Q 22 ~ Q 22 ■ Now, we examine 

expression for the excitation wave 

a similar 


E E =y 0 x + ay+vz = k xd x + k yd y + k zd z 

(B-6) 

The inverse of Equation B-2 leads to 



kyd ~ Q\l7 0 + Ql2 a + 0.32 v 



= QmY 0 + Qi2 a + Q 32 v 

(B-7) 

;d 

<7 

= 1 - (7-p v 

8 



k-zd ~ Q\2>Y 0 + Ql3 a + Q33 v 

~Q 23~r + Q33 v 

h 

(B-8) 

Now, we can identify the excitation terms in Equation B-5 and re-write it as follows 



t? 1 q a , p 2/rk , ^ , C' 2 Tuk ^ 

R kxk x ^~(ky d )y + (k zd , C?23) z 

g n 

(B-9) 

where 

(B-10) 

ra _ 3 ±a , , -± d <T-2nk ~ 

*xk ^kQ\\ + " , )Ql\ + v Q3\ 

n h 

is the response wavenumber in the x direction. If we use Equations 23 and B-7 and note that 
k yd ^ yd ' then 


o) — ky d QR + kBQ 

(B-l 1) 

and Equation 43 of the 

main text follows immediately. 
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APPENDIX C 

SYMBOLS AND SUMMARY OF COORDINATES AND WAVENUMBERS 


Symbols 


a 

c 

Cd 

c 

d 

fk 

R 

h 

i 

k 

k x ... 
m 

P 

s 

Se 

t 

W> 

W 0 

x,y,z 

B 

D 

F* 

f 

k,,k 2 ,k 3 

L 

M 

M x 

My 

Mr 

Po 

Q 

R 

AR 

T 

U 

V 0 

V 

W 

W(K,(o) 


speed of sound 

airfoil chord, measured normal to leading edge, see Figure 3 
airfoil chord, measured at constant radius, see Figure 4 

= c/R 

stagger distance, see Figure 3 
see Equation 6 

blade gap measured at constant radius 
normal distance between airfoils, see Figure 3 
imaginary unit 
scattering index 

see end of this appendix for summary of various wavenumber coordinate systems 
circumferential mode order 
acoustic pressure 

blade gap in Glegg’s coordinates, see Figure 3 

see Equation 4 

time 

upwash velocity component, see Equation 1 
complex amplitude of upwash, see Equation 1 

see end of this appendix for summary of various spatial coordinate systems 
blade count in rotor or vane count in stator 

Glegg’s potential function, described in Section 2, derived in ref. 5 
see Equation 90 

acoustic intensity component normal to cascade face, + upstream, - downstream 

y 0 R, aR, vR . Also, k/R, k 2 R, k 3 R 

normalized turbulence integral scale, A/R 

chordwise Mach number of mean flow, U/a 

Mach number in axial direction 

Mach number in tangential direction (stator frame) 

rotor rotational Mach number, QRJa 

Ambient pressure 

matrix of transformation from duct coordinates to Glegg’s coordinates 

an equivalent radius used for scaling, effective radius for noise generation 

radial extent of band of turbulence causing noise 

a random starting time used in derivation 

mean flow velocity in chordwise direction, see Figure 3 

resultant of U and W. Total mean flow velocity 

vane count 

mean flow velocity in spanwise direction, see Figure 3 
Fourier transform of gust field, see Equation 17 


I Acoustic power flux vector 

K wavenumber vector 

U mean flow vector in stator frame 
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a Glegg’s transverse wavenumber 

a ac 

a' a shifted for integration, see Equation 58 

p Vl - M 2 

y 0 Glegg’s chordwise wavenumber 

fo r» c 

(jr acoustic velocity potential, + upstream, - downstream 

$22 3D turbulence spectrum 

k see Equation 5 

K e see Equation 5 

X axial distance between chops of a turbulent eddy by a rotor 

see Equation 4 
see Equation 6 

v Glegg’s spanwise wavenumber 

v vc 

6 stagger (and flow) angle at constant radius 

# s flow angle in stator coordinates, see Figure 6 

p 0 density of mean flow 

p' acoustic density perturbation 

o interblade phase angle in Glegg’s system 

co radian frequency in stator frame 

o) g Glegg’s shifted frequency, see Equation 3 

(Ok Doppler frequency, see Equations 22 and 23 

co radian frequency in rotating frame, (same as co for stator analysis) 

co coR/a 

co coc/a 

see Equation 5 

X Glegg’s cascade stagger angle, see Figure 3 

Xe see Equation 4 

W used with subscripts for lean and sweep angles, see discussion in Section 3 

O. rotor angular speed 

A integral scale of turbulence 

<I >22 non-dimensional form of 3D turbulence spectrum, see Equation 68 

4 - 

n~ sound power upstream (+) and downstream (-), integrated over all frequency 

n* sound power spectrum upstream (+) and downstream (-), see Equation 55 
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